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Abstract 

In this article we will apply the first- and second-order supersymmetric quantum mechanics to 
obtain new exactly-solvable real potentials departing from the inverted oscillator potential. This 
system has some special properties; in particular, only very specific second-order transformations 
produce non-singular real potentials. It will be shown that these transformations turn out to be the 
so-called complex ones. Moreover, we will study the factorization method applied to the inverted 
oscillator and the algebraic structure of the new Hamiltonians. 
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1 Introduction 

Let us consider the following Hamiltonian 

h2 d2 1 2 2 
H = - — -— + -mu 2 x 2 , (1) 

2m dx 2 2 w 

where m has units of mass and w of frequency. In order to simplify, we are going to use natural units, 
such that h = m = 1, to obtain 

Moreover, by choosing appropriately the value of ui, three essentially different cases can be obtained: 

1 harmonic oscillator, 
ui = ^ free particle, (3) 
i inverted oscillator. 

These are three rare examples of exactly-solvable potentials in quantum mechanics. The first one, the 
harmonic oscillator, is a very well known system from which the technique of creation and annihilation 
operators and the whole formalism of the factorization method come from. The second, the free particle, 
has also been largely studied. This simple system allows us to work close to the limits of quantum theory, 
for example, with non-square-integrable wavefunctions with plenty of physical applications such as the 
plane waves. The third case is not so familiar: it is called either inverted oscillator, repulsive oscillator, 
inverse oscillator, or parabolic potential barrier. Although it started as an exercise from Landau's book [l], 
its physical applications have grown since the appearance of Barton's PhD thesis (published in (2|), v.g., 



* email: dbermudez@fis. cinvestav. mx 
t email: david@fis. cinvestav. mx 



1 



Factorization method and new potentials from the inverted oscillator 



Bermudez, Fernandez C. 



as an instability model, as a mapping of the 2D string theory j3j, or as a toy model to study early time 
evolution in inflationary models [4j. 

It is interesting to observe that both oscillator potentials, harmonic and inverted, are simultaneously 
produced inside an ideal Penning trap, typically used to confine charged particles [5j|6]- In its standard 
setup, a quadrupolar electrostatic field creates a harmonic oscillator potential along the symmetry axis of 
the trap, inducing confinement along this direction. In addition, in the orthogonal plane a two-dimensional 
inverted oscillator arises, driving the particles towards the trap walls. In order to compensate for the last 
effect, a static homogeneous magnetic field along the symmetry axis of the trap is also applied, but for 
zero magnetic field the two kinds of oscillator potentials are created inside the cavity. 

Mathematically, the harmonic and inverted oscillators are very much alike, and we will show that 
the solutions of one can be obtained almost directly from the other one; nevertheless, we should remark 
that physically these two systems are very different. For example the harmonic oscillator has a discrete 
non-degenerate equidistant energy spectrum with square-integrable eigenfunctions, while the inverted 
oscillator has a continuous spectrum varying from — oo to +00, which is double degenerate, and whose 
eigenfunctions are not square-integrable. 

On the other hand, a standard technique for generating new exactly-solvable potentials departing 
from a given initial one is the supersymmetric quantum mechanics (SUSY QM) (for recent reviews 
see [7| |12| ). Its simplest version, which makes use of differential intertwining operators of first-order, has 
been employed for generating Hamiltonians whose spectra differ from the initial one in the ground state 
energy level. In addition, the higher-order variants, which involve differential intertwining operators of 



orders larger than one 13 17 , allow as well the modification of one or several excited state levels. 



The SUSY techniques of first- and higher-order have been successfully applied to the harmonic oscil- 



lator 18 19 and the free particle 20 21 for generating plenty of exactly-solvable potentials. However, as 
far as we know, neither the first- nor the higher-order SUSY QM have been employed taking as a point 
of departure the inverted oscillator. In this paper we aim to fill the gap by applying the supersymmetric 
transformations to the inverted oscillator. In order to do that, in Section 2 we will get the general solution 
of the stationary Schrodinger equation (SSE) for the Hamiltonian ^ with an arbitrary energy E, which 
will remain valid even for £gC. In addition, the solutions which have a physical interpretation for the 
inverted oscillator will be identified. In Section 3 we are going to explore the factorization method for 
both systems, obtaining the bound states for the harmonic oscillator and also several sets of mathematical 



polynomial solutions, a class of solutions which have been of interest along the time (see e.g. 22-24]). 
In Section 4 we will work out the first-order SUSY QM for the inverted oscillator, while in Section 5 
we will apply the second-order one in three different situations: real, confluent and complex cases. The 
last one will be the most important case of the paper, as it is the only one that works out to obtain 
new real non-singular potentials. In Section 6 we will explore the algebraic structure for the non-singular 
potentials generated through SUSY QM and their associated eigenfunctions. Finally, in Section 7 we will 
present our conclusions. The Appendix contains the derivation of the orthogonality and completeness 
relationships for the set of eigenfunctions of the inverted oscillator Hamiltonian (see also [25]). 

2 General solution of the stationary Schrodinger equation 

First of all, let us solve the stationary Schrodinger equation for the Hamiltonian ^ with an arbitrary 
real energy E, although it is still valid for E E C: 

/ 1 d 2 uj 2 x 2 \ 

Htp(x) =(--— + —J ^(x) = Ei/,(x). (4) 
The substitution x = cj _1 / 2 y with </>(y) = i/>(w _1 / 2 i/) leads to 



Id 2 y 2 \ ^ A _ E 
2dy 



a +V W) = 77#»). ( 5 ) 
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which is the SSE for the harmonic oscillator potential in the variable y associated with the 'energy' 
E/uj [26] . Here we assume that w / 0; otherwise we immediately get the free particle problem. Thus, 
the general solution to the SSE Q reads 



ip(x) 



ux 2 /2 



E 1 



2cj 2 



-]ujx z +Dx 1 F 1 - 



E 3 



2oj 2 



, —;ujx 



(6) 



where iFi is the confluent hypergeometric function and C, D G R are constants. Note that, in general 
if>(x) $ £ 2 (IR), i.e., it does not belong to the space of square-integrable wavefunctions in one dimension 
for an arbitrary E € C. 

The analysis of the solution for the harmonic oscillator (take u) = 1 in Eq. ([6])) is widely known and 
can be found for example in 10 26 . On the other hand, since the inverted oscillator is not commonly 
studied, we will work it out next in detail. 

The general solution of the inverted oscillator is obtained from Eq. ([6]) for uj = i. For simplicity, it 
will be expressed in terms of solutions with a definite parity, the even (ip e ) and odd (ip ) solutions, which 
are given by 



^ {x)=xe~™ 2 l 2 1 F 1 



l + 2iE 1 

4 '2 ;l: 
3 + 2iE 3 



a ix 2 /2 w 



\-2iE 1 



4 ' 2' 
3-2iE 3 
4 ' 2 ; 



(7) 

(8) 



This decomposition will simplify our mathematical work in Sections 4 and 5. Therefore the general 
solution is the following linear combination 



ip(x) = Cip e (x) + Dip (x). 



(9) 



It is clear now that tp(x) is a real function, ip(x) = tp(x), for any E,C,D 6 R. 

Next, let us analyze the leading asymptotic behaviour of the solutions given in Eqs. (p7|-[8]) . To do this, 



we need the asymptotic behaviour of for \z\ 3> 1 27 

T(b) 



iiq(a, b, z) 



r(a)' 



T{b-a) 

Hence, the asymptotic behaviour for ^ e and ip Q can be straightforwardly obtained for x 1: 

n l/2 e -nE/4 



tjj e (x) 



4>o{x) 



x l/2 
7r l/2 e - OT E/4 



e i(Tr/S-x 2 /2) x -tE g- 



i(n/8-x 2 /2) x iE 



2arV2 



r(l/4-iE/2) 



i(3ir/8-x 2 /2) x -iE 

r(3/4-iE/2) 



r(l/4 + i£?/2) 

e -i(3ir/8-2; 2 /2) a ,i£; 

r(3/4 + iE/2) 



(10) 



(11) 



(12) 



By taking into account these equations, a complex linear combination is found such that now the 
terms going as x~ tE get cancelled: 



</>+(*) = N E 



2e^/ 4 T(3/A-iE/2) . 
MX) ~ r(l/4-i£?/2) ^ (a) 



(13) 



where the 'normalization' factor iVe will be chosen in order to form an orthonormal set of eigenfunctions 
in the Dirac sense (see Appendix and |25| ). On the other hand, another linearly independent solution 
can be found from ipg(x) through the reflection x — > —x: 



r(i/4-^/2) 



(14) 
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Value of E 


Classical behaviour 


Quantum behaviour 




E>0 
E = 
E < 


Goes over 

Is trapped an infinite time 
Is reflected 


Most probable goes over, some is 
Is trapped a finite time (sojourn 
Most probable is reflected, some 


reflected 
time) 
goes over 



Table 1: The classical and quantum behaviour of a particle under the inverted oscillator potential. 



It is worth noting that {ij)%{x), a = ±, — oo < E < oo} is a complete orthonormal set of eigenfunctions 
for the inverted oscillator Hamiltonian satisfying the following orthogonality and completeness relations: 



(15) 
(16) 



W£,^)= / ^ E {xW E -{x)dx = 8 a ^8{E-E'), 

J — oo 

poo 

/ dE1iZ(x$B(x') = 8(*-J), 



where 8 at<7 i and 8(y — y') denote the Kronecker and Dirac delta functions respectively (for a derivation 
of these equations we refer the reader to the Appendix and |25| ). 

One can construct now a real linear combination with a specific physical interpretation [28|29 , namely, 



e 47r / 4 r(3/4- 



iE/2) 



r(l/4 + iE/2) 



™/*T(3/4-iE/2) 
T(l/4:-iE/2) 



(17) 



The subscript L is employed due to ipi,(x) represents a particle incident from the left since its probability 
amplitude for x < is substantially bigger than the one for x > when E < 0. On the other hand, 
another linearly independent real eigcnfunction for the same E with physical meaning can be obtained 
from iPl(x) through the reflection x — > —x: 



^r(x) = IpL(-x) = 1pe(x) + 



e"/ 4 r(3/4 + iS/2) 
T(l/4 + iE/2) 



e -W4r(3/4-iE/2) 
r(l/4- iE/2) 



ip (x). 



(18) 



The subscript R denotes the fact that ipn{x) represents now a particle incident from the right. 

The classical and quantum behaviour of a particle under the inverted oscillator potential is summarized 
in Table [T] In Fig. [l] we have plotted as well the eigenfunction ?Pl(x) for different values of the energy E, 
i.e., the quantum behaviour of a particle incident from the left, which illustrates the results of Table [l] 
Since the complete set of orthonormal eigenfunctions of H {tp E (x), a = ±, — oo < E < oo} has been 
found, one can conclude that the energy spectrum of the repulsive oscillator Hamiltonian consists of the 
full real line, each eigenvalue fieK being doubly degenerated. 



3 Factorization method 

The Hamiltonian ^ is hermitian for uj — 1,0, i and it admits the well known factorization method. 
Indeed, by using the analogues of the annihilation and creation operators for the harmonic oscillator, 

^ = ^(t^+-), (19) 

it can be shown that (see also [3(5] ) 

H = a+aZ + a^a+ - |. (20) 

Since for w = 1 the operators and are mutually hermitian conjugate, (a^Y = , then for the 



harmonic oscillator the two factorizations of Eq. (20) are essentially different. On the other hand, for 
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H 




Figure 1: Eigenfunctions ipL,{x) of the inverted oscillator Hamiltonian given by Eq. (17) for the energies 
EE {-2,-1,0,1}. 



lu = i the operators af are antihermitian, (af)^ = —af, which implies that for the inverted oscillator 



the two factorizations in Eq. (20) are indeed the same (since H is hermitian). 

It is clear now that the set of operators {H, a+,a~} satisfies the following algebra 



uj, 



which leads immediately to 



Ha^j(x) = (E±oj)a^(x), 
where VK 2 -) satisfies Eq. Q, and applying n times we get 

H(at) n ^(x) = (E±m;)(atrHx). 



(21) 
(22) 

(23) 
(24) 



Note that for both cases, u) = 1 and uj = i, the general solution of Eq. Q is given by Eq. ([6| in the 
extended domain £eC; however, these solutions are not always physically admissible. Next, we will 
examine in more detail each of the two algebras. 



3.1 Harmonic oscillator algebra 



For u) = 1, Eqs. ( 19 24 ) simplify to the Heisenberg-Weyl algebra of the harmonic oscillator potential 31 
The general solution of the SSE is given by Eq. (|6| for uj = 1, from which we can obtain its bound states 
or its pure point spectrum. 



An alternative way is to take the first factorization of Eq. (20), H — a^a x + 1/2, and look for the 
extremal state ipo( x ) which is annihilated by a^, 



a 1 ip {x) = 



ip (x) = 7 r- 1 / 4 exp(-a; 2 /2). 



(25) 
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(«) 



(6) 



2 
1 

' 
-1 
-2 

-3 



•-'2 4>1 <t>0 
o o o 


$0 ^1 ^2 
• • • 


e 2 e t e 


E Ei E 2 



<p 2 <■ 




d>7~ <. 




<f>0 c 






■4 






0+c 


>4 



Re(e) 





Re(e) 



Figure 2: The complex plane for which a general solution V£ £ C can be found for both the harmonic 
and the inverted oscillators. In (a) we show the bound states (black dots) and the polynomial solutions 
(white dots) for the harmonic oscillator. In (b) we show the scattering states on the real axis (black line) 
and the polynomial solutions (white dots) for the inverted oscillator. 



In addition, ipo(x) is an eigenfunction of H with eigenvalue Eg = 1/2, and if we apply iteratively the 
creation operator we will get the remaining bound states 

( a +)n 1 
M*) = Hh^Mx) = 7^779 TJa7 m7? ™p(~x 2 /2)H n (x), (26) 



(n!)!/2 v ^ > 2 n / 2 7T 1 /4( n !)i/ 

associated with the eigenvalues E n = n + 1/2, n = 0, 1, . . . , where H n (x) are the Hermite polynomials. 
This kind of algebra is known as spectrum generating algebra. Similarly, we can obtain another extremal 
state 4>o(x) from the second factorization, H = a^a^ — 1/2, as 

a+0 o (a;)=O cf> (x) = 7r~ 1/4 exp(a; 2 /2), (27) 

which is a solution of the SSE associated with E = — 1/2 = eo. We have added the constant factor 7r~ 1//4 



by symmetry with Eq. (25), even though the wavefunction <j>a{x) is not normalizable and therefore it is 
not a bound state. We have stressed this fact by choosing a different notation </>q for this wavefunction. 

Note that the extremal state 4>o(x) can also be obtained by acting on the irregular wavefunction 
<fo{x) : defined as the second linearly independent solution of the SSE for Eq = 1/2. Although this 



wavefunction is not normalizable, it has some physical applications 32 . 

Now, we will apply iteratively the annihilation operator in order to obtain a new ladder of wavefunc- 
tions expressed in terms of Hermite polynomials, although with a different argument as compared with 
the standard case. Hence, the ladder of nonphysical wavefunctions is given by 

(a~) n i~ n 

= PjT/?^) = 2 V2 7r l/4 (w!) l/2 eMx 2 /mn(lx), (28) 

associated with the discrete nonphysical energies e„ = — n— 1/2, n = 0, 1, . . . , which nevertheless provides 
an additional set of polynomial solutions for the oscillator potential. The factor i~ n appears naturally 
when we factorize the Hermite polynomial with the correct argument. The position of the eigenvalues 



for the bound states ( 26 ) and the nonphysical energies associated with the polynomial solutions ( 28 1 on 
the complex _E-plane can be seen in Fig. [2^ a) . 

3.2 Inverted oscillator 



For u! — i Eqs. (19 24) define a zeroth order complex deformation of the Heisenberg-Weyl algebra. The 
general solution of the corresponding SSE is given by Eq. ^ with oj = i in the full complex E'-plane. 
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In this case there is no solution which fulfills the bound state condition, as it is expected because the 
inverted oscillator potential has neither local nor global wells which would support such bound states. 
Nevertheless, it is expected to have dispersive states that, under certain conditions, can be given a physical 
interpretation. Note that the most familiar system with these characteristic states is the free particle. 

As it was shown in Section 2, the natural dispersive states of the inverted oscillator are those given 
by Eqs. (13 -14) and (17 18) for any Eel. However, we can ask ourselves further: Are there solutions 



related with the Hermite polynomials for the inverted oscillator? In order to answer, let us find now the 
extremal states associated with the factorizations of Eq. (20) for ui = i. 
Let us consider in the first place 



>o(*)=0 



7r- 1 / 4 exp(-ia; 2 /2), 



(29) 



for e = i/2. Furthermore, if we apply iteratively the creation operator af we will get the following set 
of polynomial solutions of the SSE 



4>n 0) 



(n!) 1 / 



A 0) 



i/2 



2™/2 7r i/4( n !)i/: 



■&cp{-ix 2 /2)H n {i 1/2 x), 



(30) 



associated with e n = i(n + 1/2), n = 0, 1, . . . 

Similarly, the use of the second factorization leads to 



a,+0+(a:) = 



(31) 



associated with £„ 



-i/2. In addition, if we apply several times the annihilation operator a i we will 



get another set of polynomial solutions 



(n!)Va 



i/2 



2 n/a 7r i/4( n i)i/: 



■cxp{ix 2 /2)H n (i 3/2 x), 



(32) 



associated with e+ = — i(n + 1/2), n = 0, 1, 

Let us stress once again that 4>ni x ) an d <fit( x ) are solutions of an SSE for complex E- values and thus 
they do not have any physical interpretation at all, neither as bound nor as dispersive states. Nevertheless, 
it turns out that they are related with polynomial solutions of the SSE. A diagram marking the positions 
in which this kind of solutions appear for the inverted oscillator on the complex 2?-plane is shown in 
Fig. lb). 

Note that there is a direct extension of this factorization method, which is closely related to SUSY 
QM. In this generalization, instead of factorizing the Hamiltonian H in two different ways in terms of a 



pair of first-order operators, as in Eq. (20), one looks for the most general first-order operators producing 



just a single factorization [33] . Thus, when the ordering of the operators factorizing H is reversed in 
general one arrives to a different Hamiltonian. This fact has been used to generate new exactly-solvable 



Hamiltonians departing from a given initial one 33 - 35 



Next, we are going to apply the first- and second-order SUSY QM to the inverted oscillator. This 
supplies us with the building bricks for implementing the higher-order transformations, since it is known 
nowadays that any non-singular transformation of order higher than two can always be factorized in 
terms of non-singular first- and second-order SUSY transformations. Note that, although this fact was 
conjectured for the first time in 14 15 , however it was proved in a rigorous way just recently 36 37 



Then, the higher-order SUSY partners of the inverted oscillator can be obtained through iterations of 
the non-singular transformations which will be discussed here. 



4 First-order SUSY QM 

In the first-order SUSY QM one starts from a given solvable Hamiltonian Hq and tries to factorize it 
directly in terms of two first-order operators. An alternative method is to depart from an equivalent 
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intertwining relationship and to solve the resulting system of equations. We will use the last procedure 
starting from the inverted oscillator Hamiltonian 



H 



1 d 2 



+ V (x) 



1 d 2 



1 



(33) 



2 da; 2 2 ' 

and looking for a first-order differential operator which intertwines Hq with a new Hamiltonian Hi 

(34) 



in the way 33 ■ 35 
where 



H x At = A+H , 



1 

75 



d 

da; 



ot\{x) 



(35) 
(36) 



In this formalism ai(x) is called the superpotential and V\(x) denotes the new potential; until now, both 
functions are still to be determined. 

Using the standard results of the first-order SUSY QM, it turns out that a\{x) must satisfy the Riccati 
equation 

a' 1 + a 2 1 =2(V -e). (37) 



Through the substitution ol\ = u'/u, Eq. (37) is transformed into 



— -u"(x) + Vo(x)u(x) = eu(x). 



(38) 



This is the SSE for the initial potential with E = e whose general solution is given by Eq. ^. Thus, it 
is straightforward to obtain the formula for the new potential as 



V 1 (x) = V (x)-a' 1 (x) = V (x) 



u'(x) 
u(x) 



(39) 



which is real for any e, C, D 6 R. It is clear now that if the transformation function u(x) has zeros, 
then the new potential V\(x) will have singularities at those points. Since the seed solution u(x) given in 
Eq. Q, associated with an arbitrary factorization energy eeR, always has zeros because it has oscillatory 
terms that cannot get cancelled (see Fig. [I]) , it follows that it is impossible to perform real non-singular 
first-order SUSY transformations for the inverted oscillator. Note that the singular transformations are 
excluded because they change, in general, the domain of the initial potential and, consequently, the initial 



spectral problem (compare, e.g., 38 ). 



5 Second-order SUSY QM 



The second-order SUSY QM typically starts from the following intertwining relationship 11 14 

H 2 B+ = B+H Q 

where now 

1 d 2 



Ho = -\ 



2 Ax 



2 +V 2 {x), 



da; 2 



9 ^dx + k ^ 



(40) 

(41) 
(42) 
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Value of c 


Type of transformation 


c> 


Real case 


c = 


Confluent case 


c < 


Complex case 



Table 2: The three types of second-order SUSY transformations. 



The aim is to determine V-zix), g(x) and h(x) supposing that some solutions of the SSE for the potential 
Vq(x) are known. First of all, the unknown functions V^(x) and h{x) are expressed in terms of Vq(x) and 
g{x) in the way 



V 2 (x) = V -g', 
q' q 2 

Next, it turns out that g{x) must fulfill the following non-linear second-order differential equation: 



2 



12 



+ 9 2 [g' + - 1 r -2V a 



c = 0, 



with c, d € K being two integration constants. In order to solve Eq. (45 1, let us use the ansatz [39] 

g' = -g 2 + 2 lg + 2£, 



(43) 
(44) 

(45) 
(46) 



where "f{x) and £,(x) are functions to be determined. Substituting this ansatz in Eq. (45), it is obtained 
that £ 2 = c and the following Riccati equation for j(x): 



7 '+7 2 = 2(U -e), 



(47) 



with e = (d + £)/2. Therefore, to determine the function 7 one just has to solve Eq. (47), which can be 
linearized by using 7 = v! /u, leading to an SSE: 



-u' + VqU = eu. 



(48) 



The employed seed solutions u are called once again transformation functions and the associated e fac- 
torization energies. The corresponding transformations, leading to final real potentials V%{x), depend as 
well on the sign of c. Thus, they can be classified in three different cases: real, confluent, and complex 
(see Table [2). 



5.1 The real case for c > 

In this case it turns out that — \ 



rc > 0, & 



-Vc, = (d + &)/2 G 



1,2, El ^ e 2 . The 



corresponding known solutions of the Riccati equation (47) are denoted as 71 (x), 72 (x), and each one of 
them leads to (see ansatz (46)) 



</(*) 

By subtracting them we arrive at 



-g 2 {x) + 2 ll {x)g(x) + 2(e 1 -e 2 ), 
-g 2 (x) + 2 l2 (x)g(x) + 2(e 2 - 



2(ei - e 2 ) 
71 ~ 72 



{log[TU( Ul , M2 )]}', 



(49) 
(50) 



(51) 
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where Ui, u 2 are two solutions of the initial stationary Schrodinger equation associated with e\, e 2 and 
W(ui,U2) = u\v! 2 — U1U2 is their Wronskian. The new potential V 2 (x) is given by 



V 2 (x) = V (x)-{log[W( Ul ,u 2 )]}", 



(52) 



i.e., to obtain a non-singular potential V 2 a W(u\,u 2 ) without zeros is required (recall that in the first- 
order case it was directly the transformation function u the one that should not have nodes). A W(u\, u 2 ) 
without zeros could be achieved if u\ and u 2 would have alternate nodes. For the inverted oscillator this 
requirement is true in some finite interval of the a;-domain but not in the full real line. Therefore, we 
cannot produce real non-singular second-order SUSY transformations for the inverted oscillator. On the 
other hand, note that the zeros of u 2 are closer to alternate in all K as ei and e 2 become closer. 
This fact hints us to use the confluent second-order SUSY QM, a well worked algorithm where the two 



factorization energies converge to a single one 21 40,41 



5.2 The confluent case for c = 



In this case ei = e 2 := e € K, and thus the ansatz in Eq. (46) leads to a Bernoulli equation for g 



g'-2 lg = -g 2 . 

Let us suppose that 7 — u' /u is known. Thus, the general solution for g becomes 

g{x) = {\og[w{x)]Y, 

with 

w{x) = w + u 2 (z)dz, 



(53) 

(54) 
(55) 

J Xq 

where wq is a constant and xq is a fixed point inside the domain of Vq. In order to obtain a non-singular 
V 2 (x) we need to employ a w(x) without zeros. Since w(x) is a monotonically non-decreasing function, 
one can use a transformation function |41| such that 



or 



lim u(x) — 0, vjf 



lim u(x) = 0, 



\u(z)\ 2 dz < 00, 



|w(z)| 2 dz < 00. 



(56) 



(57) 



This means that u(x) should have a null asymptotic behaviour in one boundary and also be square- 
integrable over the corresponding semi-bounded interval. 

Now, from the asymptotic behaviour of ip e and tp given in Eqs. we can see that the leading 



terms fall off as 



1-1/2 



remarkably they do not depend on e G M); therefore they are in the frontier 



(but outside) of £ 2 (M), as any improper Dirac base. This means that they fulfill the first conditions of 
Eqs. (56) and (57) but not the second ones, i.e., there is no linear combination of solutions which can 



be made square-integrable over a semi-bounded interval. This fact implies that the w(x) of Eq.(55) will 
always have one zero on the real axis, and consequently the confluent second-order SUSY transformation 
will always be singular. Thus one can conclude that, although at first sight it seemed possible to perform 
the non-singular SUSY transformation through the confluent algorithm, however the transformation 
functions u{x) obtained from Eq. ^ by making E = e do not fulfill neither Eqs. (56) nor Eqs. (57). 



5.3 The complex case for c < 

For c < it turns out that e\,e 2 G C and t\ = e 2 := e. In this case we are only interested in those 
transformations producing a real potential V 2 (x), which implies that 71 (x) — ^ 2 {x) :— j(x). Following a 
procedure similar to the real case, we obtain [4211431 



j(x) = {\og[w(x)}}' 



w 



(58) 
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w(x) 



(59) 



where now 

W(u,u) 

Furthermore, it is easy to show that w(x) = w(x), i.e., w(x) £ K. 

Once again, w(x) must not have zeros in K to avoid singularities in V2(x). Since w'(x) = \u(x)\ 2 , then 
w(x) is a monotonically non-decreasing function. Thus, to assure that w(x) 7^ V x £ M. it is sufficient 



to fulfill either the conditions of Eqs. (56) or those of Eqs. (571, although both cannot be accomplished 



since the transformation function u{x) used to implement the SUSY algorithm cannot be physical, i.e., 
it is not square-integrable. The complex second-order transformations such that u(x) fulfills one of these 
conditions can produce real potentials which are always strictly isospectral to the original one. 

As it was discussed at the end of Section 5.2, for real factorization energies e <E E we cannot find 
a transformation function u(x) that fulfills either Eqs. (56 1 or Eqs. (57). On the other hand, from the 



asymptotic behaviour for ip e {x") and ipo{x) given in Eqs. ( 11|12 ), we can see that both functions contain 
two terms, one going as x~ 1 ^ 2 ~ le and other as x" 1 / 24 "* 6 . Then, for Im(e) ^ we will have one term 
falling faster and other slower than x^ 1 ^ 2 . Moreover, if a linear combination of the two solutions which 



cancels the slower terms can be found, thus a solution that fulfills either the conditions of Eqs. (56 1 or of 



Eqs. (57) is obtained. 



It is straightforward to obtain now the linear combinations that fulfill the condition of Eqs. (56), the 
resulting functions being given by 



2e-W4r( 3 /4-»e/2) / ( ^ 
u P (x,e) = rp e {x) r(l/4-ie/2) 

( \ 1 I \ 2eW4 r(3/4 + , e/2) 
u N (x,e) = rp e {x) r(l/4 + ze/2) 



(60) 
(61) 



where the labels P and iV refer to 'positive' and 'negative', according to the sign of Im(e) which gives the 
right behaviour, and we have expressed explicitly the dependence of those functions from the complex 
factorization energy e. In Fig. [3] it is shown a graph of |up(a;,e)| 2 as compared with cx x^ 1 . In this 
case the conditions of Eq. ([56]) are fulfilled and thus we can perform the non-singular complex SUSY 



transformation. Note that the transformation function Up(x, e) of Eq. (60), which is square-integrable in 



for Im(e) > 0, coincides with the solution found in Section 20.10 of 



Next, we need to construct the function w(x) defined by Eq. (59), which is a real-defined function. In 



44 



order to obtain a real non-singular SUSY partner potential [42], this w{x) should be built up by using only 
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Figure 4: The function w(x) for the complex 'energy' e = 5 + i (green line). The zero of this function lies 
in +00. 

either Up(x,e) or uj^(x,e) as transformation functions, which makes that the non-decreasing monotonic 
function w(x) vanishes either at x — > 00 if Eqs. ( [56] ) are obeyed or at x — !• —00 if Eqs. (57) are satisfied. 
In Fig. [i] we show the function w(x) built from up(x,e) for a specific factorization energy. 

This analysis shows that we can implement successfully the non-singular complex SUSY transforma- 
tions through any of the two seed solutions up(x, e) or u/v(x, e). Nevertheless, we should remember that 
each of these functions works only for half of the complex plane e, i.e., the imaginary part of e should be 
positive in order to use up(x,e) or negative for Un(x, e). However, by noticing that 



up(x,e) = u N (x,e), 



(62) 



and looking at the algorithm more carefully it is seen that for Im(e) > both the complex factorization 
energy e (with its corresponding transformation function up(x,e)) and its complex conjugate e (with 
up(x,e) = Uff(x,e)) are really used by the transformation. This means that all possible non-singular 
transformations are already covered by using either up(x,e) or u^(x, e) with e lying on the right domain, 
as it is shown in Fig. [5] 

Note that the real line of the complex plane e is excluded from the domain of non-singular SUSY 
transformations since our previous analysis is valid just for Im(e) 7^ 0. Furthermore, the points e = 
±i(m + 1/2), m = 0,1,2,... are also excluded because at those points the solutions reduce to the 
polynomials of Section 3.2 and the asymptotic behaviour either of ip e or ip a , given by Eqs. pdfl2| ), is no 
longer valid. 

The analytic expression for the new potential V2 can be found by substituting Eqs. ( 58|59 ) into 
Eq. ( 43 ) in order to obtain 

V 2 (x) =- X 4-( - { ™f) . (63) 

In Fig. [| we show several supersymmetric partners of the inverted oscillator potential, built through 
this complex algorithm. We want to remark that, in these cases, there are no new bound states at all 
being created by the transformations. 



6 Algebra of the new Hamiltonians 

In the last Subsection we have finally obtained new potentials through SUSY QM departing from the 
inverted oscillator. In fact, they constitute a two-parametric family of potentials, whose parameters are 
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TT 

Up[x, e) < 
< 


T 
1 

> 

• 

► 


< 

u N (x,e) , 

III 


E 

> 

? 

• 

> 

IV 



-2-10 1 2 



Re(e) 

Figure 5: Domain of the transformation functions Up[x, e) and un{x, e) in the complex plane e. The real 
line (Im(e) = 0) and the points e = ±i(m + 1/2), m = 0, 1,2, .. . are excluded. For quadrants I and II 
we should use up(x, e) as transformation function, and un{x, e) for quadrants III and IV. 




Figure 6: Inverted oscillator potential (blue dashed line) and its second-order SUSY partners generated by 
using Up{x, e) as transformation function for different complex values of e, {10~ 5 + 5i, (l + i)/5, 1CP 2 +i}, 
which correspond to the green, magenta, yellow lines respectively. 
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Ho 



Hn 




MipE+i(x) E + i 



■ i '/•;(•'') E 



Figure 7: Diagram of the complex second-order SUSY transformation. Hq has ladder operators a { , and 
H 2 has also ladder operators denoted by Lf , which are built from and af . 



Re(e), Im(e). We have also shown that e and e induce the same transformation and, consequently, the 
same potential V 2 (x). 

Remember that H 2 is isospectral to Hq for the only second-order SUSY transformation that works, 
namely, for the complex case. Note that B~ = (B + y , then it turns out that 

B+B- = (H 2 -e)(H 2 -e), (64) 
B~B+ = (H -e)(H -e). (65) 

Let us define now a pair of new operators Lf as 

Lf=B + afB-, (66) 

which act onto the eigenfunctions of the new Hamiltonian H 2 as ladder operators because they satisfy 
the following algebra 



[H 2 ,Lf} = ±iLt, 

[Li,L+] = P 5 (H 2 +i) - P 5 (H 2 ) = Q 4 (H 2 ), 



(67) 
(68) 



where Qa,Pb are complex polynomials of orders 4 and 5 respectively. They are given by the analogue of 
the number operator for H 2 : 



L+L- = (H 2 - e)(H 2 -e)(H 2 -e- i)(H 2 -e-i)(H 2 - i/2) = P 5 (H 2 ) 



(69) 



Note that the ladder operators Lf of H 2 are also antihermitian, (Lf)^ = —Lf, as it happens for the 
inverted oscillator. 

Next we can obtain the analytic expression for the eigenfunctions of the new Hamiltonian H 2 . Using 



Eq. (40) we will get 10 



B+tP e = y/(E-e)(E-e)ip 
B-^ } = y/(E-t)(E-e)il>i 



(2) 
E ' 



(70) 
(71) 



where iPe{%) d enotes an eigenfunction o f H associated with a n a rbitrary real energy E (it can be the 
ipf(x) of Eqs. ( fl3p4| , the ip L (x) of Eq. the ip R (x) of Eq. Q or a linear combination of both) and 
iPe\ x ) denotes the corresponding eigenfunction of H 2 (see diagram in Fig. [7]). Indeed, by substituting 
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"i — * 




Figure 8: SUSY partner potential (solid curve) of the inverted oscillator (dashed line) and four of its 
eigenfunctions associated with the energies E G {—2,-1,0, 1}. The factorization energy involved in this 
transformation is e = 10 -5 + 5i. 



in the definition of B + (see Eq. (42 1) the expressions for g(x) and h(x) of Eqs. ( 44|58 1 and after several 
simplifications we get 



w'(x) 
w(x) 



u'(x) 
u(x) 



+ 2(e-E)1> E (x). 



(72) 



We should recall that u(x) is a complex transformation function, e € C is the corresponding factorization 

(2) 

energy and E £E K is the energy associated with the initial and transformed eigenfunctions ^(a;), tp E ' (x). 
Note that the asymptotic behaviour for tp^\x) is the same as for the initial eigenfunction. In Fig. |s]a 
new potential V% (x) and four of its associated eigenfunctions are shown. 



7 Conclusions 

In this work, the supersymmetric quantum mechanics has been successfully employed to generate new 
exactly-solvable potentials departing from the inverted oscillator. We have shown that the first-order as 
well as the real and confluent second-order SUSY transformations always produce singular potentials. 
On the other hand, through the complex second-order SUSY QM it is possible to generate new real 
non-singular exactly-solvable potentials. 

It has been shown that this transformation can be achieved through two specific complex seed solu- 
tions which are the only ones fulfilling the appropriate conditions in order to produce new non-singular 
real potentials isospectral to the inverted oscillator. Furthermore, we have obtained a simple analytic 
expression for the eigenfunctions associated with the new Hamiltonian Hi- Let us note that the potentials 
generated in this paper can be used as models in every physical situation where the inverted oscillator has 
been employed before (see [2}|4"p0"] ) . This is because the new Hamiltonians are isospectral to the original 
one and the form of their potentials is quite similar. These facts, in particular, could be important to 
foresee alternative models for describing the small imperfections appearing when a real Penning trap is 
built up [5 45 46 , specially if they do not change the spectrum of the corresponding ideal arrangement. 
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Furthermore, we have studied the general algebraic structure of the original system with arbitrary uj, 
which is reduced to the harmonic oscillator case for u = 1. We have analyzed in detail as well the case 
with u) = i, which is related to the inverted oscillator and turns out to have a Heisenberg-Weyl deformed 
complex algebraic structure. For the new Hamiltonians obtained through the complex second-order SUSY 
QM applied to the inverted oscillator we have examined as well the related algebra in some detail. 

In the future we would like to analyze further the algebras associated with a more general system 
with an arbitrary u>, different from 1 and i, since the symmetry in the wavefunctions tp n , 4> n , (f>~ , 

) suggests the existence of a common structure that, perhaps, would allow us to 
understand deeper this family of oscillators. 
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of Eqs. (26|28|30|32 
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Appendix 



In this Appendix we are going to derive the orthogonality and completeness relations 

/oo 
^W E '(x)dx = 6(E - E')8 a , a ,, 
-OO 

p OO 

/ dE1%(x$E(x') = 8(x-J), 



(73) 
(74) 



for the eigenfunctions of the inverted oscillator Hamiltonian given in Eqs. Jl3- 14 1, where 

^(1/2-^/4 giB/2-1 Ti i/ 2 _ iE) 



N F = 



ir 1 / 2 r(3/4-iE/2) 



Sa- ttT i is the Kronecker delta function in the indices a and a' and 8{y — y') is the Dirac delta function in 
the index y. We will point out just the main steps of the derivation given by Wolf ( [25], Section 8.2). 

First of all, departing from the standard Fourier transform, we introduce the bilateral Mellin transform 
f^ M (A) of the function f(x) and its inverse by means of 



/OO 
OO 
p OO 



iX-l/2 



CT = ±, 



where 



x for x > 0, 
for x < 0, 

for x > 0, 
—x for x < 0. 



By substituting now Eq. ( 75 1 in Eq. ( 76 ) the following orthogonality relation is obtained: 



(27T)- 1 / dx x- iX - 1 ' 2 x l a x r 1/2 - 8{\ - X' )8 a , a , 



(75) 
(76) 

(77) 
(78) 

(79) 



On the other hand, if Eq. (76) is substituted in Eq. (75 1 one arrives at the following completeness 
relation: 



/OO 
dXxf-^x'- a iX - 1/2 = 6(x-x'), 



(80) 



meaning that the set of functions {(27r) -1 / 2 x iX , cr = ±, — oo < A < oo} constitutes a generalized 
(Dirac) orthonormal basis for £ 2 (M). 

Basically, what is left to do now is to express the eigenfunctions ip%(x) of the inverted oscillator 
Hamiltonian in terms of this orthonormal basis, so it is much easier to prove Eqs. (173 74 1. We can 



accomplish this through a clever use of the Fourier transform (see Section 7.5.12 of [25]): 

poo 

r E (x) = 2 iB / 2 (2tt)- 1 / dpp- 1 ^ 2 e ^l^ W ' 2 \ 



(81) 



where p a , a — ± are defined similarly as x a in Eqs. (77 -78 1 



Finally, the substitution of Eq. (81) inside the integral of Eq. (73) and the use of Eq. (791 leads to 



the right hand side of Eq. ( 73 ) and hence to the orthogonality relation. In the same way, by substituting 



Eq. ( 81 ) into the left hand side of Eq. ( 74 ) and using Eq. ( 80 1 one arrives to the right hand side of Eq. ( 74 ) 



and thus to the completeness relation. 



□ 
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